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Abstract 

We propose a new, renormalizable approach to nucleon-nucleon scattering in chiral effec- 
tive field theory based on the manifestly Lorentz invariant form of the effective Lagrangian 
without employing the heavy-baryon expansion. For the pion-less case and for the formula- 
tion based on perturbative pions, the new approach reproduces the known results obtained 
by Kaplan, Savage and Wise. Contrary to the standard formulation utilizing the heavy- 
baryon expansion, the nonperturbatively resummed one-pion exchange potential can be 
renormalized by absorbing all ultraviolet divergences into the leading S-wave contact in- 
1 teractions. We explain in detail the differences to the heavy-baryon formulation and present 

numerical results for two-nucleon phase shifts at leading order in the low-momentum ex- 

(N ■ 

. pansion. 

o 



1 Introduction 



The last two decades have witnessed a renewed interest in the nuclear force prob- 
lem and nuclear physics thanks to the development and application of effective field 
theory (EFT) methods. Much of this research has been influenced by the seminal 
work of Weinberg [Q3| who was the first to apply chiral perturbation theory (ChPT) 
to nucleon-nucleon scattering. Using the heavy baryon (HB) formulation of ChPT, 
he showed that reducible time-ordered nucleon-nucleon (NN) diagrams yield en- 
hanced contributions to the scattering amplitude as compared to naive dimensional 
analysis. The enhancement can be traced back to the appearance of pinch singu- 
larities emerging from the two-nucleon intermediate states. The enhanced contri- 
butions can be most easily and efficiently resummed by solving the Lippmann- 
Schwinger (LS) equation. The description of low-energy nucleon dynamics, there- 
fore, naturally reduces to the conventional, quantum mechanical A-body problem 
where the nuclear forces are defined as a kernel of the corresponding dynamical 
equation and can be derived order-by-order in ChPT. 
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Starting from the pioneering work of Ref. [2], this approach has developed rapidly 
over the last decades and is nowadays widely employed in studies of low-energy 
few- and many-nucleon dynamics and nuclear structure calculations, see H3I4I5I for 
recent review articles. While offering many attractive features such as simplicity 
and the ability to use well-developed machinery to treat few- and many-body dy- 
namics, Weinberg's approach suffers from being rather intransparent with regard to 
renormalization. One issue is related to the fact that iterations of the truncated NN 
potential within the LS equation generate contributions to the amplitude beyond the 
order one is working. These higher-order terms generally involve ultraviolet (UV) 
divergencies which cannot be absorbed by counter terms (contact interactions) in- 
cluded in the truncated potential so that one needs to employ a finite UV cutoff A 
of the order of a natural hard scale, say A ~ A x ~ M p j6l. While subleading and 
higher-order corrections to the potential do not have to be resummed in Weinberg's 
power counting scheme and can be treated perturbatively, the LS equation for the 
leading-order (LO) potential already turns out to be not renormalizable (in the usual 
sense). In particular, infinitely many counter terms are needed to absorb UV diver- 
gences emerging from iterations of the LO long-range potential due to one-pion 
exchange (OPE), whose singular l/r 3 -piece generates UV divergencies in all spin- 
triplet partial waves. This problem, in fact, shows up in every spin-triplet partial 
wave. To be specific, consider the lowest-order potential in Weinberg's approach, 

= ~4Fg 1 ' 2 q 2 + Ml +Cs + CtG1 ' ° 2 ' ( ) 

where gA, F n and M n are the nucleon axial-vector coupling, pion decay constant 
and the pion mass, respectively, <?, (Tj) denote the spin (isospin) Pauli matrices 
of the z-th nucleon and q = p 1 — p is the nucleon momentum transfer. It is easy 
to see by dimensional arguments that the 2n-th iteration of this potential in the LS 
equation generates a logarithmic divergence oc (Qm) 2n 0, where m is the nucleon 
mass and Q denotes the generic soft scale corresponding to external three-momenta 
of the nucleons and M w . In the 1 S channel, where the singular tensor part of the 
OPE potential vanishes, the coefficients in front of the logarithmic divergences do 
not involve external momenta and can be absorbed by derivative-less contact op- 
erators with multiple insertions of M 2 . The potentially enhanced contributions of 
these higher-order A/ 2 -dependent operators might become an issue if one is inter- 
ested in the quark mass dependence of nucleon-nucleon scattering but do not affect 
the predictive power of the theory in terms of describing the energy dependence 
of the phase shift at the physical values of the quark masses. On the other hand, 
in spin-triplet channels, the coefficients of the logarithmic divergences do involve 
powers of external momenta!]] Their removal requires the inclusion of an infinite 
number of higher-derivative contact interactions. For example, calculating the dia- 
grams of Fig. Q] in dimensional regularization with n spatial dimensions one finds 



1 If cutoff regularization is employed, which is normally the case in calculations with non- 
perturbative pions, one also needs to keep track of power-like divergences. The coefficients 
in front of power-like divergences generally also involve powers of external momenta. 
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(a) (6) 



Fig. 1. Two examples of iterations of the lowest-order Lippmann-Schwinger equation. 

the divergent parts proportional to m 2 M%/(n — 3) for graph (a) and, among many 
other divergent terms, m 6 (p 6 + p' 6 ) /{n — 3) for graph (b) Q. The appearance of 
such divergences seems to indicate that the same enhancement, which is responsible 
for non-perturbativeness of the OPE potentiaOl also applies to higher-order short- 
range operators. This feature is sometimes referred to as inconsistency of Wein- 
berg's approach. It should, however, be understood that this issue does not affect 
the predictive power of Weinberg's approach with regard to describing the energy 
dependence of the scattering amplitude provided a suitably chosen finite cutoff is 
employed along the lines of Ref. [6J. The predictive power of such a framework can 
be understood in terms of the modified effective range expansion H8I9B and relies on 
the knowledge of the long-range tale of the interaction driven by the pion exchange. 
We also emphasize that there is no consensus on the relevance of the inconsistency 
issue for renormalized contributions to the scattering amplitude, see e.g. Ref. ifTOll . 

The possibility to remove the UV cutoff A from the LS equation by enforcing the 
limit A — y oo (or, equivalently, A ^> A x ) non-perturbatively was also explored 
by several authors, see e.g. IfTTil . It is possible to obtain a finite, manifestly non- 
perturbative solution of the LS equation with a singular l/r 3 -potential by including 
one/no contact operator in each attractive/repulsive channel IfTTI . Unless all UV 
divergences emerging from iterations of the LS equation are absorbed by counter 
terms, such a procedure is expected to violate the low-energy theorems and is in- 
compatible with the principles of EFT |fT2|. Presently, there is still no consensus 
in the community on the most consistent and efficient way to organize the chi- 
ral expansion in the few-nucleon sector, see H6I11I12I13I14I15I16I17I18I19I20H for 
samples of different views and formulations. 

Treating the exchange of pions perturbatively as suggested by Kaplan, Savage and 
Wise (KSW) [21] obviously allows one to avoid the above-mentioned inconsistency 
and, at the same time, provides a transparent and analytically solvable EFT frame- 
work for nucleon-nucleon scattering. In the KSW approach, the OPE potential is 
shifted to next-to-leading order (NLO). The LO NN amplitude, therefore, becomes 
renormalizable perturbatively as well as non-perturbatively. Corrections beyond LO 
are treated as perturbations. While the KSW scheme is free of any inconsistencies 
with respect to renormalization, it was shown not to converge in low spin-triplet 

2 One way to justify the need to resum the OPE and the LO contact operators is to treat 
the nucleon mass as a separate scale in the problem, which is assumed to be much larger 
than the chiral symmetry breaking scale HI. 
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partial waves H22I23I24L see however Ref. 11251 for a modified formulation. 



In this paper we identify the origin of non-renormalizability of the LO NN am- 
plitude in Weinberg's approach with the non-relativistic expansion of the nucleon 
propagators. It is by now well established in the single-baryon sector of ChPT that 
the chiral power counting can be maintained without relying on the non-relativistic 
or HB expansion, see Ref. ||26l and references therein. We propose here a formu- 
lation of chiral EFT for NN scattering with non-perturbative pions based on the 
manifestly Lorentz-invariant form of the effective Lagrangian which is consistent 
and renormalizable. We will refer to this formulation as the modified Weinberg ap- 
proach. We demonstrate in the case of perturbative pions how renormalization is 
carried out in this framework and present cutoff-independent results for NN phase 
shifts at LO in the modified Weinberg approach. Our paper is organized as follows: 
In section |2] we describe the framework and specify the dynamical equation for the 
NN amplitude. Based on the (modified) KSW approach at NLO in the EFT expan- 
sion, we demonstrate in section [3] the consistency of this scheme with respect to 
the power counting. The results for non-perturbative pions at LO in the modified 
Weinberg approach are presented in section|4l Finally, our findings are summarized 
in section |5l 



2 The framework 



To derive the dynamical equation for NN scattering in chiral EFT we follow closely 
the procedure of Ref. [27] (but refrain from using high-derivative regularization and 
performing expansion in inverse powers of the nucleon mass m). We start with the 
manifestly Lorentz-invariant effective Lagrangian for pions and nucleons. It is orga- 
nized in a derivative and quark-mass expansion and consists of the purely mesonic, 
pion-nucleon and nucleon-nucleon parts, whose lowest-order contributions can be 
found e.g. in Refs. H27I28I29I , see also |[2"6ll and references therein. Following 
Weinberg 113011 , we employ time-ordered perturbation theory without performing 
non-relativistic expansion to calculate NN scattering amplitude. We decompose the 
numerator of the standard fermion propagator as 

lj) + m = 2 m P + + {j> — m ifi ) , (2) 

where P + = (1 + if>)/2 with v = (1,0,0,0), and identify the second terms as 
a higher-order correction (to be included perturbatively). For NN scattering the 
two-nucleon intermediate states generate enhanced contributions (T). Therefore, 
defining the two-nucleon-irreducible diagrams as the effective potential V, the NN 
off-mass- shell scattering amplitude T satisfies the integral equation written sym- 
bolically as 

T = V + V G T , (3) 
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where G is the two-nucleon propagator. Substituting the expansions of V, G and 
T in a small parameter (pion mass or small momenta) in Eq. ©, the leading-order 
contribution emerges from solving the equation non-perturbatively 

T = Vo + V G T , (4) 

while corrections are calculated perturbatively using the solution to the LO equa- 
tion. 

The physical amplitude is obtained from the off-shell amplitude T via U3U4 T u\u 2 
where Z$ is the residue of the nucleon propagator and Wj, Ui are Dirac spinors. To 
determine the physical amplitude order-by-order we expand the Dirac spinors in 
small quantities as 

u = u + Ui + u 2 + ■ ■ ■ , u = u + Ui + u 2 + ■ ■ ■ , (5) 

where 

«o(p) = p + u (p) , «o(p) = «(p) p + , ' ' • ■ ( 6 ) 
Consequently, the lowest-order on-shell amplitude requires the knowledge of the 
quantity T = P+P+ T P+P+ which fulfills the integral equation 

T = V + V Gf 0l (7) 

where V = P + P + V P + P + is the projected potential. In Weinberg's approach, the 
projected lowest-order potential consists of derivative-less contact interactions and 
the OPE piece. Following Ref. [27 J, we choose to treat corrections to the static 
OPE potential perturbatively so that the LO potential takes the form of Eq. (OQ) 
after switching to the two-component Pauli spinors. In the center-of-mass frame 
with incoming (outgoing) three-momentum of the nucleons p (p '), the LO equation 
takes the form 



T (p',p) = V (p',p) - 




where we have suppressed the spin and isospin indices. For the half-off- shell kine- 
matics, po is given by p = y/p' 2 + m 2 . This equation is, of course, not new and 
was proposed for the first time in Ref. BT1 . This and many other equations of a sim- 
ilar type, all satisfying relativistic elastic unitarity, have been extensively studied in 
the literature in the context of three-dimensional reductions of the Bethe-Salpeter 
equation, see e.g. Ref. ||32l . 

It is obvious that Eq. © has a milder UV behavior than the LS equation, which 
emerges from keeping the leading term in the 1 /m-expansion of the integrand. It 
is perturbatively renormalizable in the sense that all UV divergences generated by 
its iterations are absorbable into the low-energy constants (LECs) Cs.t in Eq. ©• 
Nonperturbatively, the UV behavior of this equation will be addressed in section 
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HI While all UV divergences emerging from iterations of Eq. ([8]) can be removed 
by the Cg^-counter terms, the resulting A-independent expression for the ampli- 
tude still violates the power counting due to the explicit appearance of the nucleon 
mass in the integrand and requires additional, finite renormalization. This can be 
achieved by subtractions which remove the positive powers of m and restore the 
chiral power counting, see Ref. 11331 for an extensive discussion. Such additional, 
finite subtractions in the case at hand only affect the lowest-order contact interac- 
tions of the Cs,t type. An explicit example of renormalization will be given in the 
next section for the case of perturbative pions. 

At low energies, the two-nucleon scattering amplitude calculated based on the man- 
ifestly Lorentz-invariant effective Lagrangian can, of course, be expanded in in- 
verse powers of the nucleon mass. In the HB approach, the 1 /m-expansion is car- 
ried out already at the level of the Lagrangian which implies that one first expands 
the integrand of Eq. §8§ and then calculates the scattering amplitude. This is le- 
gitimate in perturbative calculations (apart from some well-known cases such as 
e.g. the scalar form factor of the nucleon, in which the HB approach fails to re- 
produce the analytic structure of the relativistic loop diagrams for soft momenta). 
We will show in the next section that the two ways of doing 1 / m-expansion lead 
to the same result for the case of perturbative pions in the KSW framework. How- 
ever, care is required when resumming the pion exchange potential in the Weinberg 
approach. The leading nonrelativistic approximation of the two-nucleon propagator 



2 (k 2 + m 2 ) (p - \/k 2 + m 2 + ie^j p 2 - k 2 + ie 

obviously has a more singular UV behavior compared to the one appearing in 
Eq. ®. As a consequence, in contrast to the approach presented above and cor- 
responding to Eq. ([8]), an infinite set of counter terms needs to be included in the 
nonrelativistic framework if one aims at removing the UV cutoff by taking the limit 

A -> oo. 



3 Perturbative pions at next-to-leading order 



It is instructive to apply Eq. ([8]) to the case of perturbative pions where all calcula- 
tions can be carried out analytically. In order to facilitate the comparison with the 
original work by Kaplan, Savage and Wise based on the HB effective Lagrangian, 
we adopt here the normalization of the amplitude of Ref. [|2T1 . The expansion for 
the S-wave scattering amplitude in the KSW approach has the form 

A = A-i + Aq + Al-\ , (10) 
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A. 







M. 




a7> 



Fig. 2. The leading and subleading contributions to the NN scattering amplitude in the 
KSW approach. The solid dots denote the lowest-order contact operators and the leading 
pion-nucleon vertex oc qa while the filled squares refer to the subleading contact terms 
proportional to p 2 or M 2 . 

where the subscript indicates the power of the soft scale Q. The leading-order con- 
tribution A-i emerges from resummation of the LO contact interactions as shown 
in Fig. [21 Using the two-nucleon Green function from Eq. ©, the LO amplitude 
has the form 



A-x 



-C n 



(ID 



l-<7 J(p)' 

where the dimensionally regularized (DR) integral I(p) is given in n spatial dimen 
sions by 



. m 2 u 3 n 

Hp) = -ir 



d n k 



2 (2tt)« J [k 2 + m 2 } \po -Vk 2 + m 2 + i 0" 
1 



8tt 2 



m 



m 



y/ m 2 + p 2 



A + 2-21n— \m z 

nm + 2iirp — 2p sinh" 1 



0(n-3), (12) 



with the divergent quantity A defined as A = — \/{n — 3)— 7 — ln(47r) and p, being 
the scale parameter of DR. Further, the (bare) LEC C is simply the properly nor- 
malized linear combination of Cs,t- Here and in what follows, we use the notation 
p = \p |, k = \k\. Renormalization of A-i is achieved by subtracting the loop 
integral at p 2 = — v 2 with v chosen to be of order 0(Q), 



Ir(p, v ) = Hp) - Hi v ) = 

and replacing C by C^(v) which yields 

A-i 



m{y + ip) 
4ir 



+ 0(p 2 , v 2 



(13) 



-Cl(y) 



(14) 



Notice that while just using DR in combination with MS would be sufficient to 
render the expressions finite, one additional finite subtraction would have to be 
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performed in order to remove from 1^ terms of order ~ m 2 (i.e. O(Q )) which 
violate the power counting. The renormalized expression for A-\ clearly agrees 
with the KSW result of Ref. 11211 modulo higher-order terms emerging from the 
1/m-expansion of Ir(p 2 , v 2 ). 

The first corrections to the scattering amplitude are generated by diagrams shown in 
the second line of Fig. [2] The renormalized contributions of the dressed subleading 
contact operators have the form 



ICfp 2 
(C R ) 2 _ ' 

(15) 

where C R = C^(v), _D R = Df(v) denote the corresponding renormalized LECs 
and the two subtraction points in Aq are set to zero for the sake of convenience. 
The details of the calculation will be given in a separate publication. While A^ 
agrees with the HB result of Ref. [21] (modulo higher-order corrections from Ir(p, u)), 
the HB result for Aq is entirely given by the second term in square brackets. Given 
the scaling of the renormalized LECs C R ~ 0(Q~ l ) and Cf ~ 0(Q~ 2 ) E0, one 
observes that the first term in the brackets is of order ~ Q 3 while the second one 
is of order ~ Q 2 . Both approaches, therefore, again lead to the same result modulo 
corrections of a higher order. 

We now discuss the contributions involving pions. The second diagram in Fig. [2] 
simply yields the S-wave projected OPE potential, 

° "4F 2 1 + 4p 2ln Ml J' (16) 

Notice that at the order we are working, there is no need to distinguish between 
the physical and the chiral-limit values of the LECs such as gA, rn and F n . The 
renormalized contribution of the third diagram reads 



Ax 



in 



m 2 [2m 2 + p 2 — 2m^/m' 2 + p' 2 J c z 



87rC R 



AV) DfM 2 
~ (Co R ) 2 



i • 



A {III) 



' Hoop 



(P) 



2FI 

9 

m 
~2 



A-i \l R {p, v) - M 2 /n oop (p) 



d n k 



1 



(27r)" [A;2 



-VWT 



le 



711 



8np 



tan 



2p 



i, M 2 + 4p 2 
+ 2 ln ^MT~ 



+ o 



[(k-p)* + M*] 
p_ M? 
m' m 



(17) 



Again, this agrees with the HB KSW result modulo terms of a higher order emerg- 
ing from 1/m-expansion of JrQo, v) and Jnoop(p)- Finally and most interestingly, 
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for the fourth diagram we obtain 



Aiv) _ 9a ,2 
~ 4F 2 1 



M*I 2loop -I R (p, vf 



(18) 



where the scalar two-loop integral has the form 



' 21oop 



m 
T 

x 



d n hd n k 2 



1 



(2vr) 2 " 



[k 2 



Po 



k\ + m 2 



Po - \/ k\ + m2 + ^ 



l(ki 



k<?) 2 



Ml] 



m 



16tt 2 



7T 



2tt 2 



(19) 



ln8 2G 7C(S) l^Mj + V , 



2p 



where G ~ 0.916 is Catalan's constant and the ellipses refer to terms of order pm, 
M w m and higher. The HB result of Ref. [2lj for this diagram in our notation has 
the following form 



HB 



9 2 A m2 
64vr 2 F 2 



Ait Ml 



In 



Ml + 4p 2 



+ z tan 



2p 
M^ 



-ipf 
(20) 



Clearly, the difference in the constant terms in the square brackets of Eqs. (1191) 
and (1201) can be compensated by a finite shift of the LEC Df. While all relevant 
terms non-polynomial in Af 2 and p 2 are exactly the same in both cases, the HB re- 
sult features a logarithmic dependence of the renormalization scale which reflects 
the overall logarithmic divergence of the integral /2ioop when the integrand is ap- 
proximated by the leading term in the 1 /m-expansion. It is, therefore, necessary 
to include the -D 2 M 2 counter term in the HB approach at the same level as the 
diagram (a) of Fig. Q3 which appears at LO in the Weinberg approach, in order to 
absorb the corresponding logarithmic divergence. This is, in fact, the first mani- 
festation of the above-mentioned inconsistency issue of the Weinberg approach. In 
contrast, the original integral /2100P is finite and fulfills the power counting with- 
out any additional subtractions. Consequently, from the renormalization point of 
view, there is no need to promote the D 2 M^-term to LO if the OPE potential is 
treated nonperturbatively at the same footing as the Co -term within the modified 
Weinberg approach proposed in this paper. The same arguments also apply to more 
complicated diagrams with higher iterations of the OPE potential. 
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Fig. 3. 3 P and 3 P 2 phase shifts at E Uh = 100 MeV versus the cutoff A obtained from 
equation © with the OPE potential. 

4 Non-perturbative pions at leading order 

We now turn to the case of nonperturbative pions and numerically solve Eq. ([8]) 
in the partial wave basis for the LO potential given in Eq. (OQ). We employ a mo- 
mentum-space cutoff A when integrating over k in order to regularize the divergent 
integrals. As discussed in the previous sections, the LO equation ([8]) is perturba- 
tively renormalizable so that one can safely remove the cutoff by taking the limit 
A — > oo in any iteration. Nonperturbatively, the UV behavior in Eq. ([8]) can be 
understood by approximating the two-nucleon propagator for k — > oo via 



1 p + y/k 2 + m 2 

(k 2 + m 2 ) (p - Vk 2 + m 2 + ie) ~ (k 2 + m 2 ) {p 2 - k 2 + ie) 



k (p 2 — k 2 + % e) 

The UV behavior of this equation in the partial-wave decomposed form coincides 
with the one of the LS equation in 2 + 1 space-time dimensions. The OPE part of the 
potential therefore behaves in coordinate space for r — >■ as ~ l/r 2 in2 + l space- 
time dimensions. It is well known that the LS equation does not possess a unique 
solution if the strength of the attractive 1/r 2 potential exceeds some critical value 
which depends on the partial wave, see ||34)| for more details. The same sort of non- 
uniqueness emerges in the context of the Skornyakov-Ter-Martirosyan equation 
H35I36B which has also been addressed from the EFT point of view [|37ll . In the 
case at hand, we found that the non-unique solutions only appear in the 3 P partial 
wave. This situation is visualized in Fig. [3] where we compare the dependence of 
the 3 P and 3 P 2 phase shifts on the UV cutoff A at the fixed energy of Pi a b = 
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Fig. 4. Phase shifts calculated at LO in the modified Weinberg approach as functions of lab- 
oratory energy in comparison with the Nijmegen ll39l (filled circles) and Virginia Tech BDl 
(open triangles) partial wave analyses. Left panel: S- and P-waves, right panel: D-waves 
and the mixing angles 

100 MeV. While the phase shift in the 3 P2 channel quickly approaches the A — > 
oo limit, the observed limit-cycle-like behavior of the 3 Po phase shift reflects the 
non- uniqueness of solution of Eq. ®. While we still let the possibility open to fix 
the solution from physical principles without the need to rely on the data, see the 
discussion in Ref. 11381 . we follow here a more pragmatic approach of Ref. Il37l . 
Specifically, we fix the solution in the 3 P partial wave by including a counter term 
of the form C 3 p pp' / A 2 and tuning the LEC C 3 p to the Nijmegen partial wave 
analysis (PWA). Notice that the residual A-dependence of C 3 po is of a logarithmic 
type at any finite order in the loop expansion. Consequently, it is easy to see by 
dimensional arguments that the iterations of this contact interaction do not require 
the inclusion of higher-order counter terms. Therefore, the removed-cutoff limit is 
indeed legitimate from the EFT point of view in this case, contrary to the situation 
when positive powers of A appear in momentum-dependent counter terms |[T2~1|. A 
more detailed analysis of this issue will be published elsewhere. 

We are now in the position to discuss results for phase shifts. We employ the exact 
isospin symmetry as appropriate at LO and use the following values for the LECs 
entering the OPE potential 

M T = 138 MeV, F n = 92.4 MeV, g A = 1.267 . (22) 

The LECs Cs, Ct and Czpq are fitted to Nijmegen 1 5'o, 3 Si and 3 P phase shifts 
at energies P] ab < 25 MeV in the limit A — > oo. The resulting, cutoff-independent 
predictions for phase shifts in S-, P- and P-waves and the mixing angles e lj2 are 
visualized in Fig.|U Given that the calculations are carried out at LO, the agreement 
with the Nijmegen PWA is rather good. The large deviation for the 1 Sq phase shift 
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^ So partial wave 


a [fm] 


v rfml 


vo Tfm 3 l 


yn [fm 5 ] 


va [fm 7 ] 


NLO KSW from Ref. US 


fit 


fit 


-3.3 


18 


-108 


LO Weinberg 


fit 


1.50 


-1.9 


8.6(8) 


-37(10) 


Nijmegen PWA 


-23.7 


2.67 


-0.5 


4.0 


-20 



Table 1 

Predictions for the coefficients in the effective range expansion of the 1 So phase shifts (low- 
energy theorems) with perturbative and non-perturbative treatment of the OPE potential in 
comparison with the values from the Nijmegen PWA (extracted using the Nijm II potential 
fl4TT42T ). 

is also observed in LO KSW and (nonrelativistic) Weinberg approach and is well- 
known to be largely cured by the inclusion of the subleading contact interaction. 
In all other channels, the deviations between the theory and Nijmegen PWA are 
consistent with the expected corrections from higher-order terms in the expansion 
of the potential and also indicate that these corrections can be taken into account 
perturbatively. 

In addition to the predicted energy dependence of the phase shifts, the proper in- 
clusion of the pion-exchange physics can be tested in theoretical predictions for the 
coefficients in the effective range expansion 

p 2l+1 cot 8i{p) = h ^rp 2 + v 2 p 4 + v 3 p 6 + v 4 p s + ... , (23) 

where a, r and Vi denote the scattering length, effective range and shape parameters, 
respectively, and / is the orbital angular momentum. The energy dependence of 
the two-particle scattering amplitude near threshold is driven by the long-range 
tale of the interaction which imposes correlations between the coefficients in the 
effective range expansion |[23l . These correlations are determined by the long-range 
interaction and may be regarded as low-energy theorems (LETs). In tables [Hand 
|2l the LETs in the KSW and Weinberg approaches are confronted with the results 
of the Nijmegen PWA for the 1 S'o and 3 Si partial waves, respectively. Since in the 
KSW approach the LO S-wave amplitude does not involve effects due to OPE, 
one needs to go to at least NLO in order to test the LETs in this framework. The 
analytic expressions for the S-wave shape parameters at NLO in the KSW scheme 
can be found in Ref. Il23l . Clearly, the modified version of the KSW approach 
discussed in section |3] yields the same results for modulo terms of order 1/m 
and higher. The LETs are known to be strongly violated in the KSW approach 11231 , 
see tables \T\ and |2] The non-perturbative treatment of the OPE potential leads to 
an improved description of the LETs in the 1 S channel. It is, however, still rather 
poor at LO which should not come as a surprise given that the long-range part 
of the OPE potential generates only a small contribution to the 1 So phase shift. 
One may, therefore, expect that the LETs are strongly affected by the two-pion 
exchange contributions in this partial wave. In the 3 Si channel, in contrast, the 
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3 S\ partial wave 


a [fm] 


t rfml 


vo Tfm 3 l 


yn rfm 5 l 


va [fm 7 ] 


NLO KSW from Ref. US 


fit 


fit 


-0.95 


4.6 


-25 


LO Weinberg 


fit 


1.60 


-0.05 


0.8(1) 


-4(1) 


Nijmegen PWA 


5.42 


1.75 


0.04 


0.67 


-4.0 



Table 2 

Predictions for the coefficients in the effective range expansion of the 3 Si phase shifts (low- 
energy theorems) with perturbative and non-perturbative treatment of the OPE potential in 
comparison with the values from the Nijmegen PWA ||4~31 . 

LETs are well reproduced at LO in the Weinberg approach. The discrepancy for 
f 2 in the 3 S\ channel should not be taken too seriously given the very small value 
of this coefficient. We further emphasize that the errors quoted for v 3i4 refer to 
the estimated uncertainty of our numerical extraction of these parameters from the 
phase shifts. 



5 Summary and conclusions 

In this paper we applied the manifestly Lorentz-invariant form of the effective 
Lagrangian to the problem of nucleon-nucleon scattering without relying on the 
heavy-baryon expansion. The LO contribution to the scattering amplitude in the 
resulting modified Weinberg approach can be obtained by solving the LS-type of 
integral equation © with the kernel given by the OPE potential and derivative- 
less contact interactions. Contrary to its nonrelativistic counterpart, this equation is 
renormalizable, i.e. all UV divergences generated by its iterations can be absorbed 
by redefinition of the two LO contact interactions. The explicit appearance of the 
nucleon mass in the propagators, however, makes it necessary to perform addi- 
tional, finite subtractions in order to restore the proper scaling of the renormalized 
contributions in accordance with the power counting. Such additional subtractions 
only affect the values of the LECs accompanying the LO contact interactions. Con- 
sequently, the LO equation is renormalizable and consistent in the EFT sense. 

In the case of perturbative pions, the new approach is shown to reproduce the 
well-known results of the HB KSW framework modulo terms of a higher order in 
the 1 / m-expansion. When pions are treated non-perturbatively as suggested in the 
Weinberg scheme, the formulation we propose, being renormalizable, offers the ap- 
pealing possibility to remove the UV cutoff in the way compatible with the princi- 
ples of EFT. We have analyzed two-nucleon scattering at LO in the modified Wein- 
berg approach. We found that the integral equation does not possess a unique solu- 
tion in the 3 Po partial wave similarly to the Skornyakov-Ter-Martirosyan equation 
for spin-doublet nucleon-deuteron scattering. One possible way to fix the solution 
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in this channel is to include the corresponding contact interaction whose strength 
is tuned to reproduce the low-energy data ll37l . The obtained cutoff-independent 
results for phase shifts at LO in the modified Weinberg scheme are in a reasonably 
good agreement with the Nijmegen PWA. The LETs for the coefficients in the ef- 
fective range expansion are shown to be fulfilled to a very good (fair) accuracy in 
the 3 Si (%) channel. 

Clearly, the LO calculations presented here should be extended to higher orders 
in the chiral expansion. Given that the main benefit of the new formulation is 
its renormalizability, it is natural to treat higher-order corrections in perturbation 
theory. Recent studies HI 81 191201 carried out within the nonrelativistic framework 
seem to suggest that a perturbative treatment of the two-pion exchange might be 
phenomenologically successful- Also the fairly small deviations between the Ni- 
jmegen PWA and the LO results, see Fig. HI seem to be consistent with the possi- 
bility of a perturbative treatment of higher-order corrections. 

The proposed formulation offers also further advantages apart from its transparency 
with regard to renormalization. In particular, one may hope to benefit from remov- 
ing the finite cutoff artifacts which are unavoidable in the nonrelativistic framework 
with non-perturbative pions. Furthermore, the avoidance of the 1 /m-expansion can 
be advantageous in situations where the momentum scale ~ y/M^m associated 
with radiative pions must be explicitly taken into account (such as e.g. pion pro- 
duction in NN collisions and the quark-mass dependence of contact interactions 
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